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We report a comprehensive investigation of a model for peeling of an adhesive tape along with a 
nonlinear time series analysis of experimental acoustic emission signals in an effort to understand the 
origin of intermittent peeling of an adhesive tape and its connection to acoustic emission. The model 
represents the acoustic energy dissipated in terms of Rayleigh dissipation functional that depends on 
the local strain rate. We show that the nature of the peel front exhibits rich spatiotemporal patterns 
ranging from smooth, rugged and stuck-peeled configurations that depend on three parameters, 
namely, the ratio of inertial time scale of the tape mass to that of the roller, the dissipation coefficient 
and the pull velocity. The stuck-peeled configurations are reminiscent of fibrillar peel front patterns 
observed in experiments. We show that while the intermittent peeling is controlled by the peel 
force function, the model acoustic energy dissipated depends on the nature of the peel front and 
its dynamical evolution. Even though the acoustic energy is a fully dynamical quantity, it can 
be quite noisy for a certain set of parameter values suggesting the deterministic origin of acoustic 
emission in experiments. To verify this suggestion, we have carried out a dynamical analysis of 
experimental acoustic emission time series for a wide range of traction velocities. Our analysis 
shows an unambiguous presence of chaotic dynamics within a subinterval of pull speeds within 
the intermittent regime. Time series analysis of the model acoustic energy signals is also found 
to be chaotic within a subinterval of pull speeds. Further, the model provides insight into several 
statistical and dynamical features of the experimental AE signals including the transition from burst 
type acoustic emission to continuous type with increasing pull velocity and the connection between 
acoustic emission and stick-slip dynamics. Finally, the model also offers an explanation for the 
recently observed feature that the duration of the slip phase can be less than that of the stick phase. 



PACS numbers: 83.60.Df, 05.45.-a, 05.45.Tp, 62.20.Mk 
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I. INTRODUCTION 

Adhesive tapes are routinely used in a variety of situa- 
tions including daily usage as stickers, in packing and 
sealing. Yet, day-to-day experiences like intermittent 
peeling of an adhesive tape and the origin of the ac- 
companying audible noise have remained ill understood. 
This may be partly attributed to the fact that adhesion 
is a highly interdisciplinary subject involving diverse but 
interrelated physical phenomena such as intermolecular 
forces of attraction at the interface, mechanics of con- 
tact, debonding and rupture, visco-plastic deformation 
and fracture jl| jand frictional dissipation which operates 
during peeling 2, 3]. Yet another reason is that most 
information on adhesion is obtained from quasistatic or 
near stationary conditions. Apart from scientific interest, 
understanding the intermittent peel or the stick-slip pro- 
cess has relevance to industrial applications as well. For 
example, optimizing production schedules that involve 
pasting or peeling of an adhesive tapes at a rapid pace in 
an assembly line requires a good understanding of stick- 
slip dynamics. Moreover, insight into the time dependent 
and dynamical aspects of adhesion is expected to be im- 
portant in design of adhesives with versatile properties 
required in variety of applications, in understanding the 
mechanisms leading to the failure of adhesive joints as 
also in understanding biologically relevant systems such 



as the gecko Q or reorientation dynamics of cells [H]- 

Adhesion tests are essentially fracture tests designed to 
study adherence of solids and generally involve normal 
pulling off and peeling. Such experiments can be per- 
formed under quasistatic or near-stationary and nonequi- 
libriuni conditions as well. The latter kind of experiments 
demonstrate the rate dependence of adhesive properties. 
It is this rate dependence and the inherent nonlinearity 
that leads to a variety of instabilities. These kinds of 
peeling experiments are comparatively easy to setup in 
a laboratory. Moreover, the set up also allows one to 
record unusually long force waveforms and AE signals 
that should be helpful in extracting useful information 
on the nonlinear features of the system. 

One type of peeling experiment that yields dynamical 
information is carried out with an adhesive tape mounted 
on a roller subjected to a constant pull velocity [1, 0]- 
Peeling experiments have also been performed under con- 
stant load conditions 0, IB- ^^'^ P^^l velocities, the 
velocity of the contact point v keeps pace with the im- 
posed velocity V. The same is true at high velocities as 
well. However, there is an intermediate regime of trac- 
tion velocities where the peeling is intermittent. Peeling 
in this regime is accompanied by a characteristic audi- 
ble noise [1, 0, Hi • It must be stressed that these two 
stable dissipative branches refer to stationary branches. 
Even so, the stick-slip dynamics observed in the inter- 
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mediate region of pull velocities has been attempted by 
assuming an unstable branch connecting the two stable 
branches. The strain energy release rate shows a power 
law for low velocities with an exponent around 0.3. The 
high velocity branch also shows a power law but with a 
much higher exponent value of about 5.5 [6]. The low ve- 
locity branch is known to arise from viscous dissipation 
and that at high velocity corresponds to fracture. These 
studies report a range of wave forms starting from saw 
tooth, sinusoidal or even irregular wave form that has 
been termed 'chaotic' [6]. More recently, the dynamics 
of the peel point has been imaged as well . 

Stick-slip processes are usually observed in systems 
subjected to a constant response where- in the force devel- 
oped in the system is measured by dynamically coupling 
the system to a measuring device. The phenomenon is ex- 
perienced routinely, for example, while writing with chalk 
piece on a black board, playing violin or walking down 
a staircase with the hand placed on the hand-rail. A 
large number of studies on stick-slip dynamics have been 
reported in systems ranging from atomic length scales, 
for instance, stick-slip observed using atomic force mi- 
croscope to geological length scales like the stick-slip 
of tectonic plates causing earthquakes [l^, [3 ■ A few 
well known laboratory scale systems are - sliding fric- 
tion Is. llOll and the Portevin-Le Chatelier (PLC) ef- 
fect [13, [ij] , a kind of plastic instabilit y o bserved during 
tensile deformation of dilute alloys to name only 

two. Most stick-slip processes are characterized by the 
system spending a large part of the time in the stuck 
state and a short time in the slip state. This feature is 
observed both in experiment and in models. See for in- 
stance [O, [13 ■ A counter example where the time 
spent in the stuck is less than that in the slip state (ob- 
served at high applied strain rates) is the PLC effect ■ 
These studies show that while the physical mechanisms 
that operate in different situations can be quite varied 3 j 
in general, stick-slip results from a competition among 
the inherent internal relaxational time scales [ij, [l^ and 
the applied time scale. In the case of peeling, one iden- 
tifiable internal relaxation time scale is the viscoelastic 
time scale of the adhesive. Other relevant time scales 
that may be operative need to be included for a proper 
description of the dynamics. All stick-slip systems are 
governed by deterministic nonlinear dynamics. 

Models that attempt to explain the dynamical features 
of stick-slip systems use the macroscopic phenomenolog- 
ical negative force response (NFR) feature as an input 
although the unstable region is not accessible. This is 
also true for models dealing with the dynamics of the ad- 
hesive tape including the present work. In this context, 
it must be stated that there is no microscopic theory that 
predicts the negative force- velocity relation in most stick- 
slip situations except in the case of the PLC effect where 
we have provided a dynamical interpretation of the neg- 
ative strain rate sensitivity of the flow stress [13, [3 (^se 
below). There are several theoretical attempts to model 
the stick-slip process observed during peeling of an adhe- 



sive tape. Maugis and Barquin , were the first to write 
down a model set of equations suitable for the experi- 
mental situation and to carry out approximate dynami- 
cal analysis. These equations were later modified and a 
dynamical analysis of these equations was reported p^ . 
However, the stick-slip oscillations were not obtained as a 
natural consequence of the equations of motion [13, IMI ■ 
Indeed, these equations are singular|22|. Subsequently, 
we devised a special algorithm to solve these differential 
algebraic equations (DAE) [l^, [2^. This algorithm al- 
lows for dynamical jumps across the two stable branches 
of the peel force function. This was followed by convert- 
ing the DAE into a set of nonlinear ordinary differential 
equations (ODE) by including the missing kinetic energy 
of the stretched tape thereby lifting the singular nature 
of the DAE [13, Hj] . Apart from supporting dynamical 
jumps, the ODE model exhibits rich dynamical features. 
However, all these studies discuss only contact point dy- 
namics while the tape has a finite width. The ODE model 
has been extended to include the spatial degrees of free- 
dom that is crucial for describing the dynamics of the 
peel front as also for understanding the origin of acoustic 
emission [2^, [2^ . 

Acoustic emission is commonly observed in an unusu- 
ally large number of systems such as seismologically rele- 
vant fracture studies of rock samples [23,[H,|2§|, marten- 
site transformation [soL [Sll . [3^ , micro- fracturing process 
331. volcanic activity [3J], collective dislocation motion 
35, [3^ etc. The general mechanism attributed to AE is 
the abrupt release of the stored potential energy although 
the underlying mechanisms triggering AE are system spe- 
cific. The nondestructive nature of the AE technique has 
been useful in tracking the microstructural changes dur- 
ing the course of deformation by monitoring the AE sig- 
nals. For instance, it is used in fracture studies of rock 
samples (28j and more recently, a similar approach has 
been used in understanding collective behavior of dislo- 
cations [3^. In both these cases, multiple transducers 
are used to locate the hypocenters through an inversion 
process of arrival times (28l . [s^. In the latter case, by 
analysing the dislocations sources generating AE signals, 
the study establishes the fractal nature of the collective 
motion of dislocations. (In contrast to these dynamical 
studies, most studies on AE [H, [U, [s^l are limited to 
compiling the statistics of the AE signals in an effort to 
find experimental realizations of self-organized critical- 
ity [37,].) However, in the case of peeling, using multiple 
transducers is far from easy and only a single transducer 
is used leading to scalar time series. In such situations, 
dynamical information is traditionally recovered using 
nonlinear time series analysis [H, [s^ H^]. However, a 
major difficulty arises in the present case due to a high 
degree of noise present and the associated difficulties in- 
volved in curing the noise content. 

Despite large number of experimental investigations 
and to a lesser extent model studies, several issues re- 
lated to intermittent peeling and the associated acoustic 
emission remain ill understood. For instance, there are 
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no models (even in the general area of stick-slip) which 
show that the duration of the stick phase can be equal 
to or even less than that of the slip phase [§] , a feature 
which is quite unlike conventional stick-slip dynamics. 
From a dynamical point of view, this is also suggestive 
of the existence of at least three time scales. The model 
represents the acoustic energy in terms of the Rayleigh 
dissipation functional that depends on the local strain 
rate of the peel front and thus is sensitive to the nature 
of the peel front dynamics. While preliminary results of 
the model [l^, [2^ based on a small domain of parameters 
were encouraging, no systematic study of the influence of 
all the relevant time scales on the dynamics of the peel 
front was carried out. In particular, while the nature of 
experimental AE signals changes with the traction veloc- 
ity, the study of the influence of pull speed on internal 
relaxational mechanisms, the consequent peel front dy- 
namics and its relationship with the acoustic energy was 
not studied either. 

The principal objective of the present study is to un- 
derstand the various contributing mechanisms to the in- 
termittent peel process and its connection to acoustic 
emission. The objective is accomplished by carrying out 
a systematic study of the influence of the three inter- 
nal relaxational time scales namely the two inertial time 
scales of the tape mass and the roller inertia, and dis- 
sipative time scale of the peel front. In particular, we 
report the influence of the experimentally relevant pull 
velocity (covering the entire range) on the peel front 
dynamics. These studies show that the model exhibits 
rich spatiotemporal peel front patterns (including the 
stuck-peeled configurations that mimic fibrillar patterns 
seen in experiments) arising due to the interplay of the 
three time scales. Consequently, varied patterns of model 
acoustic signals are seen. Another consequence of the in- 
clusion of the three time scales is that it explains the 
recent observation that the duration of the slip phase 
can be larger than that of the stick-phase Inter- 
estingly, the model studies show that it is possible to 
establish a correspondence between the various types of 
model acoustic energy profiles with certain peel front pat- 
terns. More importantly, the study shows that even as 
the acoustic energy dissipated is the spatial average of 
the local strain rate, it can be noisy suggesting the possi- 
ble deterministic origin of the experimental acoustic sig- 
nals. Here, we report a detailed analysis of the statistical 
and dynamical analysis of the experimental AE signals. 
The study shows that while the intermittent peeling is 
controlled by the peel force function, acoustic emission 
is controlled by the dynamics of the peel front patterns 
that determine the local strain rate. This coupled with a 
comparative study of a comprehensive nonlinear time se- 
ries analysis (TSA) of the experimental AE signals for a 
wide range of traction velocities supplemented by a sim- 
ilar study on the model acoustic energy time series pro- 
vides additional insights into the connection between AE 
signals and stick-slip dynamics. In particular, the model 
displays the recently observed experimental feature that 




FIG. 1: (a) (Color online) A schematic representation of the 
experimental setup, (b) Plot of (piv") as a function of v" . 

the duration of the slip phase can be more than that of 
the stick phase with increase in the pull velocity. Finally, 
the model studies together with the dynamical analysis of 
the model acoustic signal provide a dynamical explana- 
tion for the changes in the nature of the experimental AE 
signal in terms of the changes in the peel front patterns. 

II. THE MODEL 

A typical experimental set up consists of an adhesive 
tape mounted on a roller. The tape is pulled at a constant 
pull velocity using a motor. A schematic representation 
of the set up is shown in Fig. [Ija). The axis of the roller 
passes through the point O into the plane of the paper. 
The drive motor is positioned at O'. Let the distance 
between O and O' be denoted by L P is the contact 
point on the peel front PQ. Let the peeled length of the 
tape PO' be denoted by L. Several geometrical features 
can be discussed using a projection on to the plane of the 
paper. Let the angle between the tangent to the contact 
point P and PO' be denoted by 9 and the angle I POO' 
by a. Then, from the geometry of the Fig. [Ija), we get 
L cos 9 = —I sin a and L sin 9 = I cos a — R where R 
is the diameter of the roller tape. Let the local velocity 
of the peel point be denoted by v and the displacement 
(from a uniform stuck state) of the peel front by u. Then, 
the pull velocity has to satisfy 

V ^ V + u + R cos 9 a. (1) 

As the peel front has a finite width, we define the cor- 
responding quantities along the peel front coordinate y 
(i.e., along the contact line) by v{y), 9{y) and oi{y). Then 
as the entire tape width is pulled a constant velocity, the 
above constraint generalizes to 

[V ~ v[y) - u{y) ~ R a{y) cos 0(y)] = 0, (2) 

where h is the width of the tape. However, we are inter- 
ested in the deformation of the peel front of the adhesive, 
which is a soft visco-elastic material. For the purpose of 
modeling, while we shall ignore the viscoelastic nature 
of the adhesive, we recognize its low modulus, i.e., we 
assume an effective spring constant kg (along the con- 
tact line) whose value is much smaller than the spring 
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constant of the tape material kt- This also implies that 
the force along PO' equilibrates fast and therefore the 
integrand in Eq. ([2|) can be assumed to vanish for all y. 
Thus, the above equation reduces to Eq. 

The present model is an extension of the ODE model 
for the contact point dynamics [24]. The ODE model 
already contains information on the inertial time scale of 
the tape mass that allows for dynamical jumps across the 
two branches of the peel force function. The extension 
involves introducing the Rayleigh dissipation functional 
to deal with acoustic emission apart from introducing the 
spatial degrees of freedom. The equations of motion for 
the contact line dynamics are derived by writing down 
the relevant energy terms consisting of the kinetic energy, 
potential energy and the energy dissipated during the 
peel process. The total kinetic energy Uk is the sum of 
the rotational kinetic energy of the roller tape and the 
kinetic energy of the stretched part of the tape. This is 
given by 

Here, ^ is the moment of inertia per unit width of the 
roller tape and p the mass per unit width of the tape. 
The total potential energy Up consists of the contribution 
from the displacement of the peel front due to stretching 
of the peeled tape and possible inhomogeneous nature of 
the peel front. This is given by 

The peel process always involves dissipation. Indeed, the 
peel force function with the two stable branches, one cor- 
responding to low velocities and another at high velocity 
arises from two different dissipative mechanisms. Apart 
from this, there is an additional dissipation that arises 
from the rapid rupture of the peel front which in turn re- 
sults in the accelerated motion of local regions of the peel 
front. We consider this accelerated motion of the local 
slip as the source responsible for the generation of acous- 
tic signals (25| . Any rapid movement also prevents the 
system from attaining a quasistatic equilibrium which in 
turn generates dissipative forces that resist the motion 
of the slip. Such dissipative forces are modeled by the 
Rayleigh dissipation functional that depends on the gra- 
dient of the local displacement rate [4l|. Indeed, such 
a dissipative term has proved useful in explaining the 
power law statistics of the AE signals during martensitic 
transformation [sol . [sil [3^ as also in explaining certain 
AE features in fracture studies of rock sample [4z| . Then, 
the total dissipation can be written as the sum of these 
two contributions 

R^i///M»))A.%4/'i^[2|f>f%, (5) 



where f{v) physically represents the peel force function 
assumed to be derivable from a potential function $(«) 
/ f{v)dv (see Ref. [§1). We denote the second term 
in Eq. ^ by TZae which is identified with the energy 
dissipated in the form of AE. In the context to plastic 
deformation, the acoustic energy arising from the abrupt 
motion of dislocations is given by TZae oc e^(r), where 
e(r) is the local plastic strain rate Following this, 

we interpret TZae as the energy dissipated in the form of 
AE signals. Note that ^ is the local strain rate of the 
peel front. As for the first term in Eq. (O, the form of 
the peel force function we use is given by 

f{v) = 402w°-34 + 171w°-i6 + 68e(^/'^-^) - 369.65w°-^ - 2. 

(6) 

We stress here that as we are interested in the generic 
properties of the peeling process, the exact form of the 
peel force function used here is not important as long 
as major experimental features like the magnitude of the 
jump in the velocity across the two branches, the range of 
values of the measured peel force function, in particular 
the values at the maximum and minimum, are captured. 

As can be seen from Eq. ([3]), there are two time scales; 
one corresponding to the inertia of the tape mass and 
the other due to the roller inertia. In addition, there is 
a third time scale, namely the dissipative time scale in 
Eq. ([5]) (second term). Thus, there are three internal 
relaxational time scales in the model. Apart from this, 
there is also a time scale due to the pull speed. Then 
the nature of the dynamics is determined by an interplay 
among all these time scales. 

It is more convenient to deal with scaled quantities. 
Consider introducing basic length and time scales which 
will be used to rewrite all the energy terms in scaled 
form. A natural choice for a time like variable is r = a;„t 
with w„ = kt/{b p). In a similar way, we introduce a 
basic length scale defined by c? = fmax/kt, where fmax 
is the value of f{v) at Vmax on the left stable branch. 
We define scaled variables by m = Xd = X{fmax/kt), 
I = Fd, L = L^d and R = R^d. The peel force func- 
tion / can be written as (j){v'') — f{v{v^))/ fmax- Here 
ti^ = v/vcUJud and V = V/vcWud are the dimensionless 
peel and pull velocities respectively with Vc = Vmax/^ud 
representing the dimensionless critical velocity at which 
the unstable branch starts. Using this we can define a 
few relevant scaled parameters Cj = {fmax/ktf{p/£), 
ka = kglP' /{kto?), 7„ = TuUJu/ [kto^), and y = ar, where 
a is a unit length variable along the peel front. The 
parameter C/ is a measure of the relative strengths of 
the inertial time scale of the stretched tape to that of 
the roller, k[j the relative strengths of the effective elastic 
constant of the adhesive to that of the tape material and 
7„ the strength of the dissipation coefficient. Then, the 
scaled local form of Eq. ^ takes the form 

X = {V' ~ v')vc + i?" -I- {sin a) a. (7) 
In terms of the scaled variables, the scaled kinetic en- 
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ergy and scaled potential energy Up can be respec- 
tively written as 



K 



ui 



1 

2C/ Jo 



+ +CfX\r) dm 



X 



dr. 



(9) 



The total dissipation in the scaled form is 



W =R 
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(t>{v''{r))dv' 



In ( dXjr) 
2 \ dr 



d f dC 



dC 
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dr \dX{r)J dX{r) dX{r) 
Using this, we get the equations of motion as 

VcV 



CfR' 
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X 



-X + ko 



d^X 



7» 

{1 + jj sin a) 



smaa 



)]/«c. 



values are from the left branch of with a small spa- 

tial inhomogeneity in X such that they satisfy Eq. ([7]) 
approximately. The system is evolved till a steady state 
is reached before the data is accumulated. 

The nature of the dynamics depends on the pull veloc- 
ity y, the dissipation coefficient 7„ and C/. We have 
carried out detailed studies of the dynamics of the model 
over a wide range of values of these parameters keeping 
other parameters fixed at i?^ — 0.35, P = 3.5, fco = 0.1 
{kt = 1000 N/m) and iV = 50 (in units of the grid size). 
Larger system size N = 100 is used whenever necessary. 



dr. 



m 

The first term on the right hand side is the frictional 
dissipation arising from the peel force function. The 
scaled peel force function, (f>{v^), can be obtained by us- 
ing the scaled velocities in Eq. ©. The nature of 4>{v'^) 
is shown in Fig. Hl^b). Note that the maximum occurs at 

= 1. We shall refer to the left branch AB as the 'stuck 
state' and the high velocity branch CD as the 'peeled 
state'. The second term on the right hand side denotes 
the scaled form of the acoustic energy dissipated. 

The Lagrange equations of motion in terms of the gen- 
eralized coordinates a{r), a{r), X{r) and X{r) are 



, (11) 

: 0. (12) 



(13) 



+ 7«f$.(14) 



However, Eqs. [Ti|) should satisfy the constraint 
Eq. ([7]). This consistency can be imposed by using the 
theory of mechanical systems with constraints p3| . This 
leads to an equation for the acceleration variable i}'*(r) 
obtained by differentiating Eq. ([7|) and using Eqs. (fT4|) . 



= [-X + ^ia^cosa-Rn^C-^r] 



(15) 



These Eqs. ^ [T3)) and ([15]) constitute a set of non- 
linear partial differential equations that determine the 
dynamics of the peel front. They have been solved by 
discretizing the peel front on a grid of N points and us- 
ing an adaptive step size stiff differential equations solver 
(MATLAB package). We have used open boundary con- 
ditions appropriate for the problem. The initial condi- 
tions were drawn from the stuck configuration (i.e., the 



III. TIME SERIES ANALYSIS OF 
EXPERIMENTAL AE SIGNALS 

One of the objectives is to carry out statistical and 
nonlinear time series analysis of experimental AE sig- 
nals associated with the jerky peel process with a view 
to understand the results on the basis of model stud- 
ies. Acoustic emission data files were obtained from peel 
experiments under constant traction velocity conditions 
that cover a wide range of values from 0.2 to 7.6 cm/s 
[43|. Signals were recorded at the standard audio sam- 
pling frequency of 44.1 kHz (having 6 kHz band width) 
using a high quality microphone. They were digitized 
and stored as 16 bit signals in raw binary files. There 
are 38 data files each containing approximately 1.2 x 10^ 
points. The AE signals are noisy as in most experiments 
on AE. 

Two characteristic features of low dimensional chaos 
are the existence of a strange attractor with self similar 
properties quantified by a fractal dimension (or equiva- 
lently the correlation dimension) and sensitivity to initial 
conditions quantified by the existence of a positive Lya- 
punov exponent. Given the equations of motion, these 
quantities can be directly calculated. However, when a 
scalar time series is suspected to be a projection from a 
higher dimensional dynamics, they are traditionally an- 
alyzed by using embedding methods that attempt to re- 
cover the underlying dynamics. The basic idea is to un- 
fold the dynamics through a phase space reconstruction 
of the attractor by embedding the time series in a higher 
dimensional space using a suitable time delay m El. 
Consider a scalar time series measured in units of sam- 
pling time At defined by [x{k), k = 1,2,3, N]. Then, 
we can construct d— dimensional vectors defined by = 
[x{k),x{k + T),---,x{k+{d-l)T)]; k^l,---,[N-{d- 
l)r]. The delay time r suitable for the purpose is either 
obtained from the autocorrelation function or from mu- 
tual information [soj . Once the reconstructed attractor 
is obtained, the existence of converged values of corre- 
lation dimension and a positive exponent is taken to be 
a signature of the underlying chaotic dynamics. In real 
systems, most experimental signals contain noise which 
in this case is high. There are several methods designed 
to cure the noise component [H, HO, IH, HE El ■ Usually, 
the cured data sets are then subjected to further analysis. 
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The correlation integral defined as the fraction of pairs 
of points and whose distance is less than r, is given 

by 



(16) 



where &{■■■) is the step function and Np the number of 
vector pairs summed. A window is imposed to exclude 
temporally correlated points [s^ . The method provides 
equivalence between the reconstructed attractor and the 
original attractor. It has been shown that a proper equiv- 
alence is possible if the time series is noise free and long 
[4^. For a self similar attractor C(r) ^ r^, where v is the 
correlation dimension [3^. Then, as d is increased, one 
expects to find a convergence of the slope dlnC {r) / dlnr 
to a finite value in the limit of small r. However, in prac- 
tice, the scaling regime is found at intermediate length 
scales due to the presence of noise. 

The existence of a positive Lyapunov exponent is con- 
sidered as an unambiguous quantifier of chaotic dynam- 
ics. However, the presence of superposed noise compo- 
nent, which in the present case is high, poses problems. 
In principal the noise component can be cured and then 
the Lyapunov exponent calculated [sol. liol l46i. l4l lisj. 
Here, we use an algorithm that does not require pre- 
processing of the data; it is designed to average out the 
influence of superposed noise. The algorithm, which is 
an extension of Eckmann's algorithm, has been shown 
to work well for reasonably high levels of noise in model 
systems as well as for short time series. The method has 
been used to analyze exp erimental time series as well (for 
details, see Ref. |50l [5l|). 

In the conventional Eckmann's algorithm [s^, a se- 
quence of tangent matrices are constructed that connect 
the initial small difference vector — to evolved dif- 
ference vectors ^i+k — ^j+k, where k is the propagation 
time. In the algorithm, the number of neighbors used is 
small typically min[2d, d+4] contained in a spherical shell 
of size e^. A simple modification of this is to use those 
neighbors falling between an inner and outer radii and 
eo respectively. Then, the inner shell is expected to act 
as a noise filter. However, so few neighbors will not be 
adequate to average out the noise component superposed 
on the signal. Thus, the modification we effect is to al- 
low more number of neighbors so that the noise statistics 
is sampled properly. (See for details (53, HIl.) As the 
sum of the exponents should be negative for a dissipa- 
tive system, we impose this as a constraint. In addition, 
we also demand the existence of stable positive and zero 
exponents (a necessary requirement for continuous time 
systems like AE) over a finite range of shell sizes eg- As a 
cross check, we have also calculated the correlation inte- 
gral and Lyapunov spectrum using the TISEAN package 
as well 391 . 



IV. DYNAMICS OF THE PEEL FRONT 

A systematic study of the dynamics of the model is es- 
sential to understand the influence of the various param- 
eters on the spatiotemporal dynamics of the peel front, 
its connection to intermittent peeling and to the accom- 
panying acoustic emission. From Eq. (jlOp . it is clear 
that the acoustic energy Rae is the spatial average of the 
local strain rate. As the peel front patterns determine 
the nature of acoustic energy, a detailed study of the de- 
pendence of the patterns on the relevant parameters and 
on the pull velocity should help us to get insight into AE 
generation process during peeling. 



A. General considerations on time scales and 
parameter values 

We begin by making some general observation about 
the various parameters and their influences. The dynam- 
ics of the model is sensitive to the three time scales (re- 
duced from four due to scaling) determined by the pa- 
rameters Cf, 7„ and V*. C/ is related to the ratio of in- 
ertial time of the tape mass to that of roller inertia (see 
below). The dissipation parameter 7„ reflects the rate 
at which the local strain rate relaxes. The pull velocity 
y determines the duration over which all the internal 
relaxations are allowed to occur. The range of C/ is de- 
termined by the allowed values of the tape mass m and 
the roller inertia /. Following our earlier studies, we vary 
/ from 10-5 to 10-2 and m from 0.001 to 0.1. Thus, Cf 
can be varied over a few orders of magnitude keeping one 
of them fixed. For model calculations, the dissipation 
parameter is varied from 0.001 to 1. (However, an order 
of magnitude estimate shows that 7„ << 1, see below.) 
The range of of interest is determined by the insta- 
bility domain which is from 1 to '-^ 12 as shown in Fig. 

mb). 

To appreciate the influence of inertial time scale of tape 
mass parameterized by C/, consider the low mass limit 
of the ODE model [2J| which has been shown to lead 
to the DAE model equations [25|. In this limit, the ve- 
locity jumps across the two branches of the peel force 
function are abrupt with infinite acceleration. However, 
finite tape mass introduces an additional time scale that 
leads to jumps in to occur over a finite time scale which 
in turn the magnitude of the velocity jumps. Indeed, the 
phase space trajectory need not jump to the high veloc- 
ity branch of 4>{v^), as we shall see. This can be better 
appreciated by considering the ODE model (that ignores 
the spatial degrees of freedom). Consider the relevant 
ODE model equations [13 (in unsealed form). 
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where a is shown in Fig. HJa) and u the displacement of 
the contact point, m is mass of the tape and k the spring 
constant of the tape. From, Eqs. ( flTl fT8|) . two inertial 
time scales can be identified, one corresponding to the 
roller inertia = {Rf /^Y^^ and another to that of the 
tape mass r2„ — {k/mY^^. (Note that k in Eq. (fT8|) of 
the ODE model corresponds to kt in the present model.) 
Thus, Cf in present model is directly related to the ratio 



of these two inertial time scales, 
we get 



Differentiating Eq. ([T]), 



R a cos 6 ~ Ra9 sin6. 



(19) 



Eq. ((T8|) is the force balance equation. In the limit 
m — > 0, we have the algebraic constraint f{v) — F{t){l — 
cos9{t)). Differentiating this equation shows that v di- 
verges at points of maximum and minimum of the peel 
force function f{v). This demonstrates that in the low 
mass limit, the orbits jump to the high velocity branch 
abruptly. Now consider Eq. (fTO|) that relates the ac- 
celeration of the peel point (v), acceleration of the dis- 
placement u i.e., il and a. This again is basically a force 
balance equation as can be seen by multiplying the equa- 
tion by the tape mass m. As the right hand side is small, 
any increase in one of these acceleration variables im- 
plies a decrease in the other variables. As low mass limit 
implies infinite acceleration of the peel front (v) across 
the peel force function, finite mass implies the velocity 
jumps across the peel force function is reduced. It is 
worthwhile to note that the effect of inertial time scale 
causing jumps across the unstable branch to occur at a 
finite time scale is a general feature. This has been rec- 
ognized and demonstrated experimentally in the context 
of the PLC effect 

Now consider estimating the order of magnitude of the 
dissipative time scale. The unsealed dissipation param- 
eter r„ is related to the fluid shear viscosity r] [4l| and 
thus an order of magnitude estimate can be obtained. 
Typical values of 77 for adhesives at low shear rates is 
1000 - 10000 Pa.s. As stress is directly related to 
shear viscosity, r„ can be estimated using typical dimen- 
sions of the peel front. It has been shown that deformed 
peel front dimension is about 100 fim 54, 55], the thick- 
ness of the adhesive is ~ 50 fim and the width of the 
peel front ^ 20 mm (width of the tape). It is easy to 
show that r„ 10^^ — 10^^ J.s. Thus, the range of ju 
is 10"^ - lO^'' taking rj ~ 1000 Pa.s. As some of the 
numbers used are material dependent, this is just an or- 
der of magnitude estimate. For model studies, the range 
of is taken to be from 1 to 0.001. However, we will 
not discuss the results for 7„ — 0.001 as these are similar 
to 0.01. 

Within the scope of the model, the model acoustic en- 
ergy given by i?Qe(r) = - (in tl^e dis- 
cretized form) depends on the nature of the local dis- 
placement rate. Based on this relation, some general 
observations can be made on the nature of Rae and its 
dependence on the peel front dynamics. From Eq. p4p . 
high 7„ implies that the coupling between neighboring 
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FIG. 2: (Color online a, b, e, f) (a, b) Snapshots for C/ = 
7.88, — 1.48, and 7^ = 1.0. (a) a smooth peel front with 
small amplitude high frequency oscillations due to finite roller 
inertia and (b) a smooth peel front, (c) Phase plot for an 
arbitrary spatial point on the peel front. Bold line represents 
(piv"). (d) Model acoustic energy plot, (e, f) Snapshots for 
Cf = 7.88, V = 1.48, and 7„ = 0.1. (Rugged peel front and 
the onset of stuck-peeled configuration.) 
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sites is strong and hence the local dynamics at one spatial 
location has no freedom to deviate from that of its neigh- 
bor. Thus, the displacement rate at a point on the peel 
front cannot differ from that of its neighbor. For the same 
reason, low implies weak coupling between displace- 
ment rates on neighboring points on the peel front which 
therefore can differ substantially. This clearly should lead 
to significantly more inhomogeneous peel velocity profile. 
Based on the above arguments, high 7„ should lead to 
smooth peel front and consequently sharp bursts in the 
model acoustic energy Rae that occurs during jumps be- 
tween the two branches of In contrast, when 7„ 
is small, Rae should be high as also spread out in time. 
However, as the exact nature of the peel front pattern 
is sensitive to the values of C/, pull velocity V and ju, 
the nature of Rae depends on all the three time scales. 
Indeed, one should expect that the more rapidly the peel 
front patterns change with time, the noisier the model 
acoustic energy should be. This is one feature that we 
hope to compare with experimental acoustic signals. 



B. Results of the model 

We have carried out extensive studies on the nature 
of the dynamics for a wide range of values of the pa- 
rameters stated above. The peel front dynamics is ana- 
lyzed by recording the velocity-space-time patterns of the 
peel front, the phase plots in the — plane for an 
arbitrary spatial point on the peel front and the model 
acoustic energy dissipated Rae- (Unless otherwise stated, 
these plots refer to steady state dynamics after all the 
transients have died out.) Here, we present a few repre- 
sentative solutions for different sets of parameters within 
the range of interesting dynamics. Our analysis shows 
that while the nature of the dynamics results from com- 
peting influences of the three time scales, the dissipation 
parameter 7„ appears to have a significant influence on 
the spatiotemporal dynamics of the peel front. 



1. Case (i), C'f = 7.88 - high (low) tape mass, low (high) 
roller inertia 

Given a value of C/ there is a range of values of (to, /). 
In this case, the set of values are: (0.1, 10""^), (0.01, 10~^) 
and (0.001, 10~^). The dissipation coefficient is varied 
from 7„ = 1 to 0.01. For high 7„ — 1.0, only smooth 
peeling is seen independent of the magnitude of the pull 
velocity. The peel front switches between the low and 
high velocity branches of the peel force function (/'(f^). 
Plots of the smooth nature of the entire peel front are 
shown in Figs. [DJa, b) for V = 1.48. Figure [^a) shows 
the nature of the peel front when the system is on the 
low velocity branch of 4>{v''), i.e., the local velocities of 
all spatial elements follow the AB branch of 4>{v'')- The 
small amplitude synchronous high frequency oscillation 
of the entire peel front results from the roller inertia. 




FIG. 3: (Color online) (a-d) Snapshots during the peeling 
process for C/ = 7.88, V = 1.48, and 7„ = 0.01 : (a) highly 
rugged peel front when the system is on the left branch of 
4>iv"), (b) the onset of peel process, (c) stuck-peeled configu- 
ration and (d) resulting nearly uniform peeled state. 

(Compare the values of v'^ in the two figures.) A phase 
plot in the ~ v'^ plane for an arbitrary point on the 
peel front is shown in Fig. [D(c). The small amplitude os- 
cillation of the peel front shown in Fig. [Ifa) corresponds 
to the velocity oscillations in the phase plot (Fig. HJc)). 
As high Cf implies relatively low values of J, it can be 
shown (on lines similar to Ref. 22]) that the orbit sticks 
to the stationary branches (slow manifold) of the peel 
force function 0(w*) jumping between the branches only 
at the limit of stability typical of relaxation oscillations. 
The corresponding model acoustic energy Raeij) shows a 
sequence of small amplitude spikes corresponding to the 
small amplitude oscillations arising from the roller inertia 
[Fig. IHd)] followed by large bursts that occur at regular 
intervals. The bursts result from the peel front jumping 
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from the stuck to the peeled state and back. Note that 
the duration of the bursts are short compared to duration 
between them. 

However, as we decrease 7„ to 0.1 keeping V = 1.48, 
we observe rugged and stuck-peeled configurations. The 
rugged pattern is seen when the system is on the AB 
branch of Even so, on reaching the hmit of sta- 

bihty, the entire contact Hne peels nearly at the same 
time as shown in Fig. [H^e). But once it jumps to the 
high velocity branch CD of (p{v''), the peel front that has 
nearly uniform peel velocity commensurate with that of 
the right branch of becomes unstable and breaks 

up into stuck and peeled segments as shown in Fig. [Hf). 
The width of these segments increases in time with a 
concomitant decrease in the magnitude of the velocity 
jumps of peeled segments, eventually the entire peel front 
goes into a stuck state. Then, the cycle restarts with the 
peel front switching between the rugged and stuck-peeled 
(SP) states. The phase plot is similar to that for 7m = 1 
again sticking to the slow manifold. The model acoustic 
energy dissipated Rae is also similar to that for 7,j = 1 
except that the large bursts are comparatively broader as 
should be expected due to presence of stuck-peeled con- 
figurations that contribute to large changes in the local 
velocity. 

As we decrease 7„ to 0.01, the observed patterns are 
similar to those for 7„ — 0.1 but the sequence of the peel 
front patterns is different. Starting with a low velocity 
configuration that is even more rugged compared to that 
for 7„ =0.1 as shown in Fig. Elja), the peel process starts 
with a small stuck segment getting peeled [Fig. [3]Jb)]. 
There after, several stuck segments peel out leading to 
a stuck-peeled pattern as shown in Fig. [3fc), eventually, 
the entire peel front peels-out leaving a nearly uniform 
peeled state as shown in Fig. EJd) (with a velocity com- 
mensurate with the high velocity branch of (jiiv^)). This 
is again destabilized with some segments of the peel front 
getting stuck as in the case of 7„ = 0.1 (similar to Fig. 
[2Jf)). The number of such stuck segments increases with 
time, eventually the whole peel front goes into a stuck 
state. The cycle restarts. The phase plot is similar to 
"fu — 0.1 and 1. Indeed, for a given C/, independent 
of 7„ the phase plot changes only when is increased. 
But Rae shows broader bursts compared to 7„ = 0.1 as 
the corresponding stuck-peeled configurations last longer. 
Even so, the duration of the SP configurations in a cycle 
is short, i.e., the duration of the bursts is short compared 
to the duration between them. 

Now we consider the influence of increasing the pull 
velocity (keeping C/ fixed at 7.88) which in turn should 
leave less time for internal relaxational mechanisms to op- 
erate. Intuitively one should expect that some patterns 
observed for low V may not be seen for higher values of 
V^. 7n = 1 case is uninteresting for the reasons stated 
above. But, reducing 7„ to 0.1 does provide some de- 
gree of freedom for the local dynamics to operate at each 
point. Even so, for = 2.48, the peel front switches be- 
tween a SP configuration with most segments momentar- 
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FIG. 4: (Color online a,e) (a, b) Snapshot of a stuck-peeled 
configuration and model acoustic energy respectively for C/ = 
7.88, ju = 0.01 and V" = 2.48. (c, d) Phase plot and model 
acoustic energy plot respectively for C/ = 7.88, V" = 4.48, 
and 7u — 0.01. (e) Snapshot of long lived stuck-peeled con- 
figuration for Cf = 7.88, V = 5.48, and 7„ = 0.01 and (f) 
the associated model acoustic energy. 
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ily in the stuck state (similar to Fig. [Sljb)) and a configu- 
ration that has several stuck-peeled segments (similar to 
Fig. EI^c)). The corresponding — phase plot shows 
that the orbit jumps moves slightly beyond the upper 
value of 4>{v'^) and jumps back from the right branch even 
before reaching the minimum of (^(w'*) ( not shown). As 
we decrease ju to 0.01, the rugged configuration seen for 
V = 1.48 is no longer seen and only SP configurations 
are observed as shown in Fig. HJa). The SP configura- 
tions are dynamic in the sense, segments that are in stuck 
state at one time become unstuck at a later time and vice 
versa. For this case (C/ = 7.88, = 2.48, 7„ = 0.01), 
these rapid changes occur over a short time scale. Con- 
sequently, the model acoustic energy is quite noisy as 
shown in Fig. [UJb) but has a noticeable periodic compo- 
nent. The points of minima correspond to configurations 
that have fewer peeled segments compared to those near 
the peak of Rae- The phase plot in the X'^ — v'^ plane 
is limited to the upper part of (j){v^ ) . Even as the phase 
plots for any two spatial points look similar, there is a 
phase difference. For instance, at any given time, the 
phase point of stuck segment will be on the left branch 
while that for peeled point will be on the right branch. 

As we increase to 4.48 (keeping C/ at 7.88), there 
is even lesser time for peel front inhomogeneities to relax 
and thus, we observe a smooth peeling for 7^ = 1 and 
as also for 0.1. As we decrease 7„ to 0.01, we see only 
SP patterns (not shown but similar to Fig. [DJa)). The 
corresponding X^ — phase plot for an arbitrary point 
on the peel front shown in Fig. [IJc) is confined to the top 
of (/>(w*). The corresponding Rae is noisy and irregular 
as shown in Fig. [UJd). However, when we increase 
to 5.48, initially, one does observe the patterns switching 
between rugged and SP configurations. If we wait long 
enough, we observe only SP configurations that are dif- 
ferent from those for lower V. In this case, the stuck 
and peeled segments are long lived. A top view of the 
SP pattern is shown in Fig. [D^e). The phase space orbit 
in the X'' — v'^ plot is pushed beyond the upper limit of 
(/)(u*). The energy dissipated is quite regular (but ape- 
riodic) unlike that for lower pull velocities as shown in 
Fig. Sl^f). This regularity is clearly due to the long lived 
nature of these SP configurations. The long lived na- 
ture of the SP configurations for high pull velocity is a 
general feature, i.e., the duration over which the stuck 
segments remain stuck (peeled segments remain peeled) 
increases as we increase the pull velocity. The dynamics 
is no longer interesting beyond V^^ = 6.48 as only smooth 
peeling is seen. 



2. Case (a), Cf — 0.788 - high (and low) tape mass, high 
(and intermediate) roller inertia 

For this value of Cf, the allowed set of values of (m, /) 
are (0.1, lO'^), (0.01, lO'^) and (0.001,10"'*). The dy- 
namics is more interesting for this case as there is a scope 
for competition among the three time scales. 



We first study the dynamics keeping V" = 1.48 and 
varying the dissipation parameter. For 7^ = 1.0, the 
uniform nature of the peel front seen for Cf — 7.88 dis- 
appears and even for short times, stuck-peeled config- 
urations are seen. The peel front patterns stabilize to 
stuck-peeled configurations as shown in Figs. [DJa, b). 
As can be seen these SP patterns have only a few stuck 
or peeled segments with moderate velocity jumps and 
smooth variation along the peel front unlike the SP con- 
figurations discussed earher. (Note that the SP configu- 
ration in Fig. [3b) has more stuck segments compared to 
Fig. ISI^a).) The moderate velocity jumps can be under- 
stood by noting that the phase space orbit never visits 
the high velocity branch of (/'(w^) as can be seen from 




time 

FIG. 5: (Color online a, b) (a, b) Snapshots of stuck-peeled 
configurations for C/ = 0.788, V" = 1.48, and 7„ = 1.0. Note 
that (a) has fewer peeled segments compared to (b). (c) Phase 
plot for an arbitrary point on the peel front, and (d) model 
acoustic energy plot. 
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FIG. 6: (Color online a, c) (a, b) Snapshot of a stuck-peeled 
configuration for C/ = 0.788, V = 1.48, and 7^ = 0.1 and 
model acoustic energy plot respectively, (c, d) Plot of 'an 
edge of peeling' configuration for C/ = 0.788, V" — 1.48, and 
7u = 0.01 and the corresponding model acoustic energy. 



Fig. [5fc). It is interesting to note that the trajectory 
stays close to the unstable branch of (f>{v'^) even after at- 
tempting to jump from the low velocity branch. Such 
orbits are reminiscent of canard type solutions [5^ . The 
trajectory is irregular and is suggestive of spatiotemporal 
chaotic nature of the peel front. The energy dissipated 
Rae shown in Fig. EUd) is continuous and irregular due 
to the dynamic SP pattern as should be expected, but 
there is a noticeable periodic component. The rough pe- 
riodicity of Rae can be traced to fact that the peel front 
configurations switch between patterns with more stuck 
segments and less stuck segments. (From the number 
shown on a;-axis, Fig. Elja) can be identified with mini- 
mum and Fig. [D^b) with the peak of Rae in Fig. [IJd). 
See the marked arrows as well.) 



As we decrease 7„ to 0.1, the SP configurations ob- 
served have more stuck and peeled segments compared 
to 7„ = 1 (compare Fig. [B]^a) with Fig. EKa)). However, 
the magnitude of the velocity jumps remains moderate 
as in the previous case. This is again due to the fact 
that the orbit never visits the high velocity branch of 
(/)(«''). (Recall that given a value of C/, the phase plot 
remains the same for different 7^ values as long as V 
is fixed). Indeed, for this value of C/ = 0.788, the or- 
bit never jumping to the high velocity branch is a conse- 
quence of finite inertia of the tape mass compared to that 
of the roller inertia as discussed earlier. For this case, the 
model acoustic energy Rae is also irregular and contin- 
uous as shown in Fig. [6jb) with a noticeable periodic 
component. Now, if we decrease ju further to 0.01, the 
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FIG. 7: (Color online a, b) (a, b) Snapshots during peel pro- 
cess for Cf = 0.788,1/" = 2.48, and 7„ = 1.0. (c) Model 
acoustic energy plot, (d) Phase plot for an arbitrary spatial 
point on the peel front. Bold line represents (j^iv"). 
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peel front pattern displays increased number of stuck and 
peeled segments with each stuck segment having only a 
few contiguous stuck points as can be seen from Fig. [6jc). 
Note also that there is a large dispersion in the magni- 
tudes of the velocity jumps of the peeled segments even 
as the largest one is significantly smaller than the value 
of CD branch of As can be seen from the Fig. EUc), 

even though the pattern is dynamic, the segments that 
are stuck are barely so. Thus, the configuration shown 
in Fig. ^c) gives the feehng of a critically poised state. 
The corresponding X* — phase plot (similar to that 
shown in Fig. [5]^c)) is irregular and possibly suggestive 
of spatiotemporal chaotic nature of the peel front. The 
acoustic energy Rae is very irregular without any trace 
of periodicity as shown in Fig. [HJd). 

We now consider the influence of increasing the pull 
velocity V^. As we increase to 2.48, the spatiotempo- 
ral patterns seen for 7^ = 1.0, 0.1 and 0.01 are slightly 
different from those for = 1.48. For 7„ = 1.0, the peel 
process goes through a cycle of configurations shown in 
Figs. [7{a, b). It is clear that Fig. ^a.) has more segments 
in the stuck state while Fig. [7]^b) is the usual kind of SP 
configuration except that the stuck and peel segments are 
fewer. For this case, the stuck and peeled segments last 
longer than those for V — 1.48. The corresponding Rae 
for each 7„ exhibits noisy bursts overriding a periodic 
component. A typical plot for 7^ = 1 is shown in Fig. 
[7jc) . From the time labels as also the arrows shown, the 
minima and maxima in Rae can be identified with Figs. 
[7]Ja, b) respectively. The orbit in the X'^ — v'' plane moves 
into regions much beyond the values allowed by 4>{v'^) as 
is clear from Fig. [7fd). The phase plots for •^u — 0.1 and 
0.01 are similar to this case. 

For 7„ — 0.1 also, the peel front pattern goes through 
a cycle of stuck-peeled configurations (with more stuck 
and peeled segments than for 7„ — 1.0) and stuck seg- 
ments (similar to Fig. [7]Ja)). Yet, the energy dissi- 
pated Rae is similar to Fig. [7Kc) for 7^ = 1.0 which 
is surprising considering that there are more stuck and 
peeled segments compared to 7„ — 1 case. This can 
be traced long lived of the stuck or peeled configura- 
tions that hardly change over a cycle (as in the case of 
Cf = 7.88, = 5.48, 7„ = 0.01, see Fig. He)). The peel 
process is similar even for 7„ = 0.01. 

As we increase the peel velocity to 4.48, the infiuence 
of this time scale on the peel front pattern is discern- 
able even for 7„ = 1.0. The spatiotemporal patterns of 
the peel front switches sequentially from nowhere stuck 
configuration shown in Fig. [5^a) to stuck-peeled configu- 
ration with few stuck and peeled segments shown in Fig. 
[8] (b). Note that there are very few stuck and peeled seg- 
ments. The corresponding Rae exhibits noisy periodic 
pattern similar to Fig. [T^c) for = 2.48. The - v' 
phase plot in Fig. [5Kc) shows that the orbit can move 
much beyond the values allowed by (^(w"). As we de- 
crease 7„ to 0.1, the nowhere stuck configuration [Fig. 
[SKa)] is replaced by a partly stuck, partly peeled config- 
uration and a SP configuration. For 7„ = 1.0 case ( Fig. 
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FIG. 8: (Color online a, b, d, e) (a, b) Snapshots during peel 
process Cf = 0.788,1/" = 4.48, and 7„ = 1.0. (c) Phase 
plot for an arbitrary spatial point on the peel front. Bold 
line represents (piv"). (d, e) Snapshots during peel process 
Cf = 0.788, V = 4.48, and 7„ = 0.01. 



IHIc)), the AT* — phase plot is slightly different as the 
orbit makes several loops before it jumps to low velocity 
branch without visiting the high velocity branch of (^(w"). 
The nature of Rae is still noisy and periodic similar to 
Fig. EIc). As we decrease 7„ to 0.01, the peel process 
goes through SP configurations shown in Figures. [5Jd, 
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e). Note that Fig. ^e) has large dispersion in the magni- 
tude of velocity jumps of the peeled segments compared 
to that in Fig. [HJd). It is worth emphasizing that the in- 
crease in the number of stuck and peeled segments with 
decrease in 7„ is a general feature. Despite the higher 
number of stuck and peeled segments, Rae for 7 — 0.01 
is similar to that for 7„ = 1.0 as these peel front config- 
urations are long lived which again is a general feature 
observed at high pull velocities. Finally, it should be 
stated that for Cf = 0.788, in general the velocity varia- 
tion along the peel front is much more smooth compared 
to other values oi Cf. The dynamics is uninteresting 
beyond = 7.48 as only smooth peeling is seen. 



3. Case (in), Cf — 0.00788, low tape mass and high roller 

inertia 

For this value of C/, there is just one set of values 
of tape mass and roller inertia, namely, m = 0.001 and 
/ = 0.01. As the tape mass is low, this also corresponds 
to the DAE type of solutions for each spatial point. Thus, 
the velocity jumps between the two branches of the peel 
force function will always be abrupt with the roller iner- 
tia playing a major role in allowing the orbits to jump 
between the branches of the peel force function as demon- 
strated earlier f22l|. 

Consider the influence of the dissipation parameter 
keeping = 1.48. For 7u = 1, peeling is uniform 
and thus the whole peel front switches between the two 
branches of the peel force function. The acoustic energy 
shows a bunch of seven double spikes that appear at regu- 
lar interval as shown in Fig. [H^a). (The number of spikes 
is correlated with the number of cascading loops seen in 
the X" — v" phase plot, see below.) As we decrease 7„ to 
0.01, the peel front goes through a cycle of patterns with 
only few peeled segments and those with large number 
of stuck-peeled segments as shown in Figs. MJ^) and (c) 
respectively. The phase plot in the X'^ — plane of an 
arbitrary point on the peel front jumps between the AB 
to CD branches of (^(w'*). As shown in Fig. [Hd), in a 
cycle, the trajectory starting at the highest value of (piv") 
stays on CD for a significantly shorter time compared to 
that on the left branch. The orbit then cascades down 
through a series of back and forth jumps between the 
two branches of (For V = 1.48, independent of 

"fu value, the nature of the phase plot is the same with 
seven loops.) 

The corresponding model acoustic energy consists of 
rapidly fluctuating time series with an overall convex en- 
velope of bursts separated by a quiescent state as shown 
in Fig. [5)Je). (Contrast this with Fig. [DJa) for 7„ — 1.) 
From the time labels in Figures. MJ^, c), both configura- 
tions belong to the region within the bursts [Fig. [5{e)]. 
To understand this complex pattern of bursts in Rae we 
have looked at the fine structure of each of these bursts 
along with the evolution of the associated configurations. 
One such plot is shown in Fig. [5](f) which shows that fine 
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FIG. 9: (Color online b, c) (a) Rae as a function for Cf = 
0.00788, V = 1.48, and 7„ = 1. (b) a stuck-peeled con- 
figuration with more stuck segments compared to (c) where 
nearly equal number of stuck and peeled segments for C/ = 
0.00788, V" = 1.48, and 7„ = 0.01. (d) The corresponding 
phase plot for an arbitrary spatial point along the peeling 
front. Bold line represents <?f>(i;°). (e) Corresponding Rae{r) 
as a functions of time and (f) substructure of (e). 
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structure consists of seven bursts within each convex en- 
velope. These seven bursts can be correlated with the 
seven loops in the phase plot shown in Fig. [Hd). The 
time interval marked LM in the phase plot corresponds 
largely to stuck configuration (not shown) and hence can 
be easily identified with the quiescent region in Rae- Fol- 
lowing the peel front patterns continuously, it is possible 
to identify the sequence of configurations that leads to 
the substructure shown in Rae [Fig. [DJf)]. For instance, 
the loop marked PQRST in the X'^ — w'* plot corresponds 
to the burst between P and T in Fig. [IJf). During this 
period, the configuration at P is largely in the stuck state 
(as in Fig. [Ifb)) which gradually evolves with more and 
more segments peeling out [Fig. M,c)] as the trajectory 
moves from P ^ Q ^ R ^ S . As the number of stuck 
and peeled segments reaches a maximum, Rae reaches the 
peek region. Then, during the interval corresponding to 
S to T the number of peeled segments decreases abruptly. 
Thereafter, the next cycle of configurations (correspond- 
ing to the next loop in the phases plot) ensues. 

As we increase the pull velocity, the peel front is 
smooth for 7„ = 1.0 as also for 0.1 for the entire range 
of pull speeds. However, for 7„ = 0.01, as we increase 
to 2.48, the peel process goes through a cycle of SP 
configurations shown in Figs. [Toya. b). Note that there 
is a large dispersion in the jump velocities as is clear 
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FIG. 10: (Color online a, b) (a, b) Snapshots during peel 
process for C/ = 0.00788, V = 2.48, and 7^ = 0.01. (c) 
Corresponding model acoustic energy plot. 




FIG. 11: (a) (Color online) Snapshot during peel process 
Cf = 0.00788, V" = 4.48, and 7„ = 0.01. (b) Model acoustic 
energy plot, (c) Phase plot for an arbitrary spatial point on 
the peel front for C/ = 0.00788, V = 4.48, and 7„ = 0.01. 
Bold line represents (jiiv") 



from Fig. [TOT a). The corresponding Rae shows rapidly 
fluctuating triangular envelope of bursts with no quies- 
cent region seen for V = 1.48 case. This is shown in 
Fig. fTOT c). The corresponding phase plot is similar to 
Fig. [DJd) but has twelve loops. In addition, the value 
of the upper loop extends far beyond that allowed by 
(^(u**). We also see a fine structure similar to that in Fig. 
[21^e). As in the previous case, it is possible to identify 
configurations that correspond to minima and near the 
maxima of Rae- As we increase further to 4.48, only 
SP configurations are seen. The energy dissipated Rae 
shows continuous bursts overriding a sawtooth form as 
shown in Fig. [TTT b'). The AT" — phase plot shows large 
excursions way beyond the peel force function values as 
shown in Fig. [TlTc). 

A general comment may be relevant regarding large 
excursions of the trajectory in the phase plot X" — as 
we increase the pull velocity. This is easily explained for 
the low Cf (low tape mass, high roller inertia). It is clear 

from Eq. ^ that m ^ we have F{t) ^ (^i+i-na(t)) ■ 
As a{t) can takes on positive and negative values, one can 
see Fmax and Fmin are determined by minimum (nega- 
tive) and maximum values of sin a as argued in |22l | . It 
is possible to extend this argument to finite tape mass 
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TABLE I: Largest Lyapunov exponent for the model for var- 
ious parameter values. For all C/, the LLE reaches a near 
value zero for — 5.48. 

case. Finally it must be stated that the dynamics is no 
longer interesting beyond V'^ — 7.48. 
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FIG. 12: Largest Lyapunov exponents of the model for Cf — 
0.788 and 7„ = 0.01: (a) V = 1.48 , (b) V = 2.48 and (c) 
V" = 4.48. 



4- Spatiotemporal Chaotic Dynamics 

As discussed above, there are several sets of parameter 
values for which the phase plots are irregular which may 
suggest the possibility of spatiotemporal chaotic dynam- 
ics. To verify this, we have calculated the largest Lya- 
punov exponent (LLE) for all the cases using the model 
equations. (The transient solutions for the first 2000 time 
units have been ignored for calculating the LLE.) Figure 
[T^ shows a plot of the largest Lyapunov exponent for 
Cf — 0.788 and 7„ — 0.01 for various pull speeds. The 
value of LLE for — 5.48 is close to zero (not displayed 
in the figure). As can be seen, the LLE is positive being 
largest for = 1.48 decreasing to near zero value for 
V = 5.48. Table U shows the values of LLE for various 



parameter values for which the spatiotemporal dynamics 
has been detected. From this we conclude that the dy- 
namics of the peel front is spatiotemporally chaotic for a 
range of parameters values. 

V. ANALYSIS OF AE SIGNALS 
A. Statistical Analysis of AE signals 

The acoustic emission data obtained from experiments 
are fluctuating and noisy only within the domain where 
the peel process is intermittent form 0.2 to 7.6 cm/s. 
As known from early experiments Q, force wave forms 
change as the traction velocity is increased. Correspond- 
ingly, the nature of the AE signals also change with the 
traction velocity. At low traction velocities, the AE sig- 
nals have a burst like character appearing at nearly reg- 
ular intervals separated by oscillatory decay of the am- 
plitudes. These bursts are correlated with the stick-slip 
events. With increasing traction velocity, the bursts be- 
come increasingly irregular and continuous. Examples of 
burst and continuous type of AE time series are shown 
in Figures. [Qfa. b). 

As shown in the previous section, the nature of the 
model AE signal depends on parameter values. In gen- 
eral Rae can be of noisy burst type overriding a periodic 
component, continuous and irregular, rapidly fluctuating 
triangular envelope of bursts or simply a set of spikes. 
Clearly, interesting cases for comparison with the exper- 
imental AE signals are those where Rae is continuous 
and noisy. Simplest quantity to compare is the nature of 
the model acoustic signal with the energy of the exper- 
imental AE signal (i.e., square the amplitude). Figures 
llSf c. d) show a comparison between model acoustic sig- 
nal for Cf = 7.88, 7„ = 0.01 and V = 1.48 and energy 
of the experimental signal for V — 0.4 cm/s. Both show 
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TABLE IL Statistical and dynamical invariants for the exper- 
imental AE signals for typical traction velocities. The second 
and third columns show power law exponents tuai and mA2 
corresponding to small and large amplitudes. When only a 
single power law is seen, mA2 is the exponent value. Fourth 
to sixth columns list the values of the correlation dimension 
V, the largest exponent Ai and Kaplan- Yorke dimension Dky 
obtained from dynamical analysis of the AE signals. NC cor- 
responds to nonchaotic dynamics where we did not find any 
convergence of the correlation dimension. 
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FIG. 13: (a, b) Raw AE signal for V = 1.6 cm/s and 7.6 
cm/s respectively, (c) Burst like model acoustic energy plot 
for Cf = 7.88, V = 1.48 and 7„ = 0.01. (d) Burst like ex- 
perimental acoustic energy for V — 0.4 cm/s. (e) Continuous 
model acoustic energy plot for Cf = 7.88, = 4.48 and 
7ti = 0.01. (f) Continuous experimental acoustic energy for 
V = 6.4 cm/s. 



burst type emission. As another example Figs. llSf e. 
f) show respectively the continuous model signal (for 
Cf = 7.88, 7„ = 0.01 and V = 4.48) and experimen- 
tal acoustic energy for V = 6.4 cm/s. 

Given an experimental time series, the simplest statis- 
tical quantity to compute is the statistics of events. The 
definition of events depends on the physical situation, 
vifhich in the case of the AE signal may be the time in- 
terval between the bursts of AE, the amplitude of bursts 
etc. Indeed the former has been computed 0. Here we 
compute the distribution of the amplitudes of the AE 
signals. The difference between the maximum and next 
minimum, denoted by AA can be taken to be a measure 
of the amplitude of the AE signal. (In experiments, it 
is measured by setting a cut-off and measuring all am- 
plitudes larger than the cut-off.) We have computed the 
distribution of the amplitudes D{AA) for all the 38 data 
files for pull velocities starting from 0.2 to 7.6 cm/s. Sur- 
prisingly, we find only power law distributions for all the 
data files, i.e., D{AA) ~ AA~™^; we do not find peaked 
distributions. For small traction velocities, we find a sin- 
gle power law crossing over to a two stage power law for 
high traction velocities. A typical single power law dis- 
tribution for V = 1.6 cm/s is shown in Fig. [Hf a) with 
an exponent niA = 2.15. Figure [HT b) shows a two stage 
power law for high velocity V = 5.0 cm/s. The exponent 



values are niA = 0.32 and 3.0 respectively for small and 
large amplitude regimes. The transition from a single to 
two stage power law distribution occurs with the devia- 
tion for small values seen in Fig. I14f b) becoming more 
dominant with increase in the pull velocity. A two stage 
power law (over one order of magnitude range) is first ob- 
served for V = 3.0 cm/s. The exponent values are func- 
tions of the pull velocity. Table [Til shows the exponent 
values for a selected set of pull velocities. Even though 
statistical features are easy to calculate, they are suffi- 
ciently discriminating. The analysis will be useful while 
comparing the cured data files as also with the statistics 
of model acoustic energy signals. 

These results may be compared with the statistics of 
the amplitude of the model energy bursts Rae, i-e., from 
the maximum to the next minimum. Denoting ARae to 
be the amplitude of Rae{T), let D{ARae) be the distri- 
bution of the amplitude of Rae- For Cf = 0.00788, only 
7u = 0.01 case is interesting (see Fig. [TUFd) correspond- 
ing to = 2.48). To determine D{ARae), we use long 
time series (typically 10^ points in units of the integra- 
tion step). For this case, we find a two stage power law 
as shown in Fig. \W[d) i. e., D{ARae) AR-^'^ with 
the exponents ttie ^ 0.60 and 2.0 for the small and large 
amplitude regimes respectively. Figure 114( c) shows an 
example of a single power law for Cf — 0.788, 7„ ~ 0.1 
and V = 1.48. The exponent value is itie — 0.7. To 
compare, we note that the experimental time series refers 
to the amplitude of the AE signals while the model signal 
Rae is the energy. Thus, the two exponents are related 
throughrTiB = (r7iA + l)/2. Using the value of shown 
in Fig. [HTd). we obtain m^i = 0.2 for the exponent cor- 
responding to small amplitude regime (of the AE signal) 
and mA2 = 3.0 for large amplitudes. Clearly, the values 
are in reasonable agreement with the exponent values for 
the experimental signals [Fig. [TWb)]. Unlike the experi- 
mental signal where scaling regime is good for all values 
of the pull velocity, for the model acoustic energy Rae, 
the distributions for ARae show a power law statistics 
(of at least one order of scaling regime) only for a certain 
sets of parameter values. Table IIIII lists the exponent 
values wherever the power law distribution is seen. For 
high and low values of C/ we find two stage power law 
distributions. However, for the intermediate C/( 0.788), 
only a single stage power law is found. It is interesting 
to note that the power law generated here is purely of 
dynamical origin. 



B. Dynamical Analysis of AE time series 

In Section IV, we showed that the peel front pat- 
terns for several sets of parameters are spatiotemporally 
chaotic. More importantly, the model acoustic energy is 
quite irregular even as it is of dynamical origin. This 
suggests the possibility that the experimental AE sig- 
nals could be chaotic. However, often time series have 
undesirable systematic component, which needs to be re- 
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FIG. 14: (a) Plot of a single stage distribution of the ampli- 
tudes of the uncured data for V = 1.6 cm/s. (b) Plot of a two 
stage power law distribution of the amplitudes of the uncured 
data for V = 5.0 cm/s. (c) A plot of a single stage power law 
distribution for the magnitudes of acoustic energy dissipated 
for Cf = 0.788, 7„ = 0.01 and V = 1.48. (d) A plot of a 
two stage power law distribution for the magnitudes of model 
acoustic energy for C/ = 0.00788, 7„ = 0.01 and V" = 2.48. 
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TABLE III: Statistical and dynamical quantities for the model 
acoustic signal. Columns 4 and 5 show the power law expo- 
nents. Columns 6 to 8 list the correlation dimension, positive 
exponent and Lyapunov dimension respectively. 



moved from the original data. For instance, in the PLC 
effect, the stress-strain time series has an overall increas- 
ing stress ari sing from the work hardening component 
of the stress j5P| which needs to be subtracted. In the 
present case, the experimental data for high pull veloci- 
ties does show a background variation. A simple way of 
eliminating this background component is to use a win- 
dow averaging and subtract this component from the raw 
data. Moreover, as stated in the introduction, the exper- 
imental AE data are quite noisy and therefore it is nec- 



essary to cure the data (using standard noise reduction 
techniques [s^) before subjecting them to further analy- 
sis. Simple visual checks for the existence of chaos such as 
phase plots, power spectrum etc. have been carried out. 
We have also used singular value decomposition, false 
neighbor search etc. Figures llSf a. b) show the raw and 
cured data respectively for V = 5.0 cm/s. Clearly, the 
dominant features of the time series are retained except 
that small amplitude fluctuations are reduced or washed 
out [1^. Statistical features like the distribution func- 
tion for the amplitude of the AE signals, power spectrum 
etc. are not altered. For instance, the two stage power 
law distribution for the amplitude of AE signals for the 
raw data (shown in Fig. [HTb)) is retained except that 
the exponent value for the small amplitude regime is re- 
duced from 0.32 to 0.24 without altering the exponent 
corresponding to large amplitudes. This reduction is un- 
derstandable as small amplitude fluctuations are affected 
during curing. 

The cured data are used to calculate the correlation 
dimension for all the data files. However, for calculating 
the Lyapunov spectrum using our algorithm, raw data is 
adequate as our algorithm is designed to process noisy 
data. (In contrast, calculating the Lyapunov spectrum 
using the TISEAN package requires the cured data). To 
optimize the computational time, all our calculations are 
carried out using one fifth of each data set as each file 
contains large number of points ^ 10^ points, and there 
are 38 data sets. 

Typical autocorrelation time is about four units in 
sampling time. However, using a smaller value of r = 1, 
we have calculated the correlation integral C(r) for all 
the data files. Converged values of correlation dimension 
are seen only in the region of pull velocities in the subin- 
terval 3.8 to 6.2 cm/s. A log-log plot of C[r) for the pull 
velocity 3.8 cm/s is shown in Fig. [TST c) for d = 9 to 13. 
A scaling regime of more than three orders of magnitude 
is seen with v ^ 2.80 ± 0.05. This is at the beginning of 
the chaotic window. 

We have calculated the Lyapunov spectrum using our 
algorithm. The Lyapunov spectrum for V = 3.8 cm/s 
is shown in Fig. [T5f dV (The outer shell radius is 
kept at 60 = 0.065.) Note that the second expo- 
nent is close to zero as should be expected of contin- 
uous flow systems. Using the spectrum, we have cal- 
culated the Kaplan- Yorke dimension (also called Lya- 
punov dimension) D^y using the relation D^y = j + 

^i;;^!' ; ELi ^» > 0;E'=i'A« < O. The value so ob- 
tained in each case should be consistent with that ob- 
tained from the correlation integral. For the case shown 
in Fig. [I5td) we get Duy = 2 + 1.70/1.81 = 2.94 consis- 
tent with V ~ 2.80. (Typical error bars on the first three 
Lyapunov exponents are ±0.01, ±0.005 and ±0.05. Thus 
the errors in Dky values are ±0.05.) 

As an example of converged value of correlation dimen- 
sion near the upper end of the chaotic domain, a log-log 
plot of C(r) for = 5.0 cm/s is shown in Fig. [TSTe) with 
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FIG. 15: (a, b) Raw and cured AE signal respectively for 
V — 5.0 cm/s. (c) Correlation integral for pull velocity 3.8 
cm/s for d = 9 to 13. Dashed lines are guide to eye. (d) 
Lyapunov spectrum for the same data file, (e) Correlation 
integral for pull velocity 5.0 cm/s from d = 7 to 10. Dashed 
lines are guide to eye. (f) Lyapunov spectrum of the AE 
signals for traction velocity 5.0 cm/s. 
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are in the range of 2.7 to 3.0. We have also calculated 
the Lyapunov spectrum using the TISEAN package using 
cured files. The Dky values obtained from the TISEAN 
package are uniformly closer to the v values, typically 
Dky = v + 0.1. Finally, we note that the positive expo- 
nent decreases toward the end of the chaotic domain (6.2 
cm/s). These results (see Table HI]) show unambiguously 
that the underlying dynamics responsible for AE during 
peeling is chaotic in a mid range of pull speeds. 

In order to compare the low dimensional chaotic nature 
of the experimental AE signals with the model acoustic 
signal, we have analyzed the low dimensional dynamics of 
Rae{T) using the embedding procedure after subtracting 
the periodic component. 

We have computed the correlation dimension and Lya- 
punov spectrum for the entire instability domain. A log- 
log plot of the C(r) is shown in Fig. [TCT a) for = 5 to 8. 
The convergence over more than three orders of magni- 
tude is clear. The value of = 2.20 ± 0.05. For this file, 
we find stable positive and zero exponents for a range of 
Co values. A plot of the spectrum for C/ = 0.00788, 7„ = 
0.01 and V = 2.48 (e,, = 0.08) is shown in Fig. 
Using this we get Dky = 2 + 0.32/0.77 2.4 which is 
again consistent with i/ = 2.2 ± 0.02. 

We have calculated both correlation dimension and 
Lyapunov spectrum of Rae for a range of values of the 
parameters. For each C/, we find converged values of v 
and Dky within a window of pull speeds. Generally, the 
range of v is between 2.15 to 2.70 while Dky is in the 
range 2.4 to 2.90. Table IIIII shows the values of corre- 
lation dimension and Dky for various sets of parameter 
values. It is interesting to note that the magnitude of the 
largest exponent for the model AE signal also decreases 
as we increase the pull velocity, a feature displayed by 
the experimental time series as well. 



VI. SUMMARY AND CONCLUSIONS 



FIG. 16: (a) Correlation integral of Rae(T) for C/ — 
0.00788, 7„ = 0.01 and V = 2.48 with = 2.2 (d = 5 to 
8). Dashed lines are guide to eye.(b) The Corresponding Lya- 
punov spectrum for Rae. 



V = 2.73 ±0.05 for d = 7 to 10. Again, the scaling regime 
is seen to be nearly three orders of magnitude. The Lya- 
punov spectrum for the data file is shown in Fig. [TsT f). 
The calculated Lyapunov dimension from the spectrum 
is Dky = 2 + 1.73/1.80 — 2.96 which is again consistent 
with 1/ = 2.73. 

The values of i' for all the files are found to be in the 
range 2.55 to 2.85 ± 0.05 as can be seen from Table HIl 
We have calculated the Lyapunov spectrum for the full 
range of traction velocities and we find (stable) positive 
and zero exponents only in the region 3.8 to 6.2 cm/s, 
consistent with the range of converged values of v as can 
be seen from Table [TTl The corresponding values of Dky 



In summary, the present investigation is an attempt 
to understand the origin of the intermittent peeling of 
an adhesive tape and its connection to acoustic emission. 
At the conceptual level, we have established a relationship 
between stick-slip dynamics and the acoustic energy, the 
latter depends on the local strain rate [4l[ which in turn 
is controlled by the roughness of the peel front. As the 
model is fully dynamical, one basic result that emerges 
is that the model acoustic energy is controlled by the na- 
ture of spatiotemporal dynamics of the peel front. Fur- 
ther, even as the model acoustic emission is a dynamical 
quantity, the nature of Rae turns out to be quite noisy de- 
pending on the possible interplay of different time scales 
in the model. Thus, the highly noisy nature of the ex- 
perimental signals need not necessarily imply stochastic 
origin of AE signal; instead, they could be of determin- 
istic origin. This motivated us to carry out a detailed 
analysis of statistical and dynamical features of the ex- 
perimental AE signals. Despite the high noise content. 
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we have been able to demonstrate the existence of finite 
correlation dimension and positive Lyapunov exponent 
for a window of pull speeds. The Kaplan- Yorke dimen- 
sion (for various traction velocities) calculated from the 
Lyapunov spectrum is consistent with the value obtained 
from the correlation integral. Thus, the analysis estab- 
lishes unambiguously the deterministic chaotic nature of 
the experimental AE signals. Interestingly, the largest 
Lyapunov exponent shows a decreasing trend toward the 
end of the chaotic window, a feature displayed by the 
model acoustic signal as well. The work also addresses 
the general problem of extracting dynamical information 
from noisy AE signals. A similar analysis of the model 
acoustic energy shows that Rae is chaotic for a range 
of parameter values. More importantly, several qualita- 
tive features of the experimental AE signals such as the 
statistics of the signals and the change from burst to con- 
tinuous type with increase in the pull velocity are also 
displayed by Rae ■ The observed two stage power law dis- 
tribution for the experimental AE signals [Fig. [HTb)] is 
reproduced by the model [Fig. El^d)]. It must be empha- 
sized that this power law distribution for the amplitudes 
is completely of dynamical origin. This result should be 
of general interest in the context of dynamical systems 
as there are very few models that generate power laws 
purely from dynamics. The only other example known 
to the authors is that of the PLC effect where the ampli- 
tude of the stress drops shows a power law distribution 
within the context of the Ananthakrishna model [l^ . 

The spatiotemporal patterns of the peel front are in- 
deed rich and depend on the interplay of the three time 
scales. Although, the nature of spatiotemporal patterns 
is quite varied, they can be classified as smooth syn- 
chronous, rugged, stuck-peeled and even nowhere stuck 
patterns. As expected on general consideration of dy- 
namics, rich patterns are observed for the case when all 
the three time scales are of similar magnitude (illustrated 
for Cf = 0.788). All spatiotemporal patterns, except 
the smooth synchronous peel front are interesting. As 
a function of time, the nature of the peel front can go 
through a specific sequence of these patterns (depending 
on the parameter values). The most interesting pattern 
is the stuck-peeled configuration which is reminiscent of 
fibrils observed in experiments (s^ . [ssl . . Even among 
the SP configurations, there are variations, for example, 
rapidly changing, long lived, edge of peeling etc. Despite 
the varied range of patterns, a few general trends of the 
infiuence of the parameters on the peel front patterns 
are worth noting. First, in general the number of stuck 
and peeled segments increases as 7„ is decreased. Sec- 
ond, as the pull velocity is increased, the rapidly varying 
stuck-peeled configurations observed at low pull veloci- 
ties become long lived. The dynamical signature of these 
two parameters are reflected in the nature of the phase 
space orbit. For instance, given a value of Cf and 7„, the 
nature of the phase space orbit changes only when V is 
increased which allows the orbit to move way beyond the 
values of (f>{v^). 



The study of the model shows that while the intermit- 
tent peeling is controlled by the peel force function, the 
dynamics of the peel front is influenced by all the three 
time scales. This together with the dynamical analysis 
of the experimental acoustic emission signals establishes 
that deterministic dynamics is responsible for AE during 
peeling. The various sequences of peel front patterns and 
their time dependences lead to quite varied model acous- 
tic signals. These can be classified as bunch of spikes, 
isolated bursts occurring at near regular intervals, con- 
tinuous bursts with an overall envelope separated by qui- 
escent state, continuous bursts overriding a near periodic 
triangular form, irregular waveform overriding a periodic 
component, and continuous irregular type. Interestingly, 
our studies show that there is a definite correspondence 
between the model acoustic energy and the nature of peel 
front patterns even though Rae{T) is the spatial average 
of the local strain rate. Despite this, two distinguishable 
time scales in Rae{T) can be detected, one correspond- 
ing to short term fluctuations and another corresponding 
to overall periodic component. The short term fluctua- 
tions can be readily identified when the model acoustic 
signal is fiuctuating without any background component 
(see for example Fig. EJd)). These rapid changes in Rae 
arise due fast dynamic changes in the SP configurations. 
The minimum in Rae corresponds to the situation where 
the average velocity jumps of the SP configurations are 
smaller compared to that at the preceding maximum. In 
contrast, the overall periodicity in i?Qe(r) (for instance 
see Fig. [TJ^c) among many other cases) can be identified 
with the changes in the peel front patterns that occur 
over a cycle in the phase plot. The minima in the Rae{T) 
corresponds to the peel patterns where more segments 
are in stuck state than in the peeled state while the max- 
ima corresponds to more stuck and peeled segments (see 
Figs. [Zta), (b)). The corresponding phase plot usually 
goes through a cycle of visits between the low and high 
velocity branches. 

Often, however, the nature of the model acoustic en- 
ergy signal can be complicated as in the case of low tape 
mass {Cf = 0.00788). Even in such cases, some insight 
is possible. This is aided by the analysis of the corre- 
sponding phase plot. For example for the low Cf case 
where the roller inertia plays an important role in the 
dynamics, the rapidly fiuctuating acoustic energy has an 
overall triangular envelope [Fig. (Hie)]. On an expanded 
scale, the convex envelope consists of seven local peaks 
[Fig. EKf)]. Each of these is generated when the various 
peel front segments make abrupt jumps between the two 
branches of the peel force function. Note that the phase 
space orbit has seven loops in this case [Fig. (H^d)]. The 
general identification of the minima in Rae with patterns 
that have more stuck segments than peeled segments still 
holds. Similarly, the maxima in Rae usually correspond 
to the presence of large number of stuck-peeled configu- 
rations. 

The above correspondence between the model acoustic 
energy and the peel front patterns provides insight into 
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the transition from burst to continuous type of AE seen 
in experiments as a similar transition from burst type to 
continuous type is also seen in the model acoustic energy 
(for large C/ = 7.88 ). At low pull velocities, the peel 
front goes through a cycle of patterns where most seg- 
ments of the peel front (or the entire peel front) spends 
substantial time in the stuck state switching to stuck- 
peeled configuration. As the duration of the SP config- 
uration is short and velocity bursts are large, Raeir) is 
of burst type [Fig. [THfc)]. With increasing pull velocity, 
only dynamic stuck-peeled configurations are seen which 
in turn leads to continuous AE signals [Figures. Ufa, 
b)]. This coupled with time series analysis of the model 
acoustic signal shows that the associated positive Lya- 
punov exponent decreases with increase in the traction 
velocity. This is precisely the trend observed for exper- 
imental signals as well. Thus, the decreasing trend of 
the largest Lyapunov exponent can be attributed to the 
peel front breaking up into large number of small seg- 
ments providing insight into stick-slip dynamics and its 
connection to the AE process. 

The present study has relevance to the general area 
of stick-slip dynamics. As mentioned earlier, models for 
stick-slip dynamics use negative force-drive rate relation. 
In such models, the phase space orbit generally sticks to 
the slow manifold (stable branches) of the force-drive rate 
function. This leads to clearly identifiable stick and slip 
phases, the former lasting much longer than the latter. 
However, recent work on imaging the peel point dynam- 
ics shows that the ratio of the stick phase to the slip 
phase, is about two or even less than unity for high peel 
velocities. While all the known models of the peel process 
predict that the duration of the stick phase is longer than 
that of the slip phase, our model displays the experimen- 
tally observed feature. This feature emerges in the model 
due to the interplay of the three time scales aided by in- 
complete relaxation of the relevant modes. Our studies 
show that only for low pull velocity and high C/ do we 
observe the stick phase lasting much longer than the slip 
phase. As the pull velocity is increased, and for all other 
parameter values, we find that the duration of the slip 
phase (peel velocity being larger than unity) is nearly the 
same as or less than that of the stick phase (peel veloc- 
ity less than unity). Further, the present model provides 
an example of the richness of spatiotemporal dynamics 
arising when more than two time scales are involved. In 
this context, we emphasize that the introduction of the 
Rayleigh dissipation functional to model the acoustic en- 
ergy is crucial for the richness of the spatiotemporal peel 
front patterns. It is important to note that this kind of 
dissipative term is specific to spatially extended systems 
as it represents relaxation of neighboring points on the 
peel front. 

The present study has relevance to time dependent is- 
sues of adhesion. For instance, apart from the fact that 
the time series analysis addresses the general problem of 
extracting dynamical information from noisy AE signals, 
it may have relevance to failure of adhesive joints and 



composites that are subject to fluctuating loads. The fail- 
ure time can be estimated by calculating the Lyapunov 
spectrum for the AE signals. If the largest Lyapunov 
exponent is positive, the inverse of the exponent should 
give an estimate of the time scale over which the failure 
can occur and hence could prove to be useful in predict- 
ing failure of joints. The present study should also help 
to optimize production schedules in peeling tapes. 

Finally, several features of the present study are com- 
mon to the PLC effect even though the underlying mech- 
anism is very different. In this case the repeated occur- 
rence of stress drops during constant strain rate defor- 
mation [3, [58 1, are associated with the formation 
and possible propagation of dislocation bands that are 
visible to the naked eye. The phenomenon occurs only 
in a window of applied strain rates. The instability is 
attributed to the pinning and unpinning of dislocations 
from solute atmosphere, yet, the dominant feature un- 
derlying the instability is the negative strain rate sen- 
sitivity of the flow stress that has two stable branches 
separated by an unstable branch. Clearly, these features 
are similar to the occurrence of the peel instability within 
a window of pull velocities and the existence of unstable 
branch in the peel force function. Further, the Anan- 
thakrishna (AK) model for the PLC instability predicts 
that the stress drops should be chaotic in a subinterval of 
the instability domain [iTj . This prediction has been ver- 
ified subsequently through the analysis of experimental 
stress-strain curves obtained from single and polycrystals 
[sol . [59I [60l . [Slj . This feature is again similar to the ex- 
istence of chaotic dynamics observed in a mid range of 
pull velocities in the peeling problem, both in experiment 
and in the model. In the case of the PLC effect, one finds 
that the positive Lyapunov exponent characterizing the 
stress-time series decreases toward the end of chaotic win- 
dow, both in experiments and in the AK model. Again 
this feature is also seen in the present peel model as also 
in experimental AE signals. Dynamically, in the case of 
the AK model for the PLC effect the decreasing trend of 
the positive Lyapunov exponent has been shown to be a 
result of a forward Hopf bifurcation (HB) followed by a 
reverse HB [l^. In the case of peeling problem as well, 
the instability begins with a forward HB followed by a 
reverse HB. Finally, in the PLC effect (both in exper- 
iments and the AK model), as in the peeling problem, 
the duration of the slip phase can be longer than that of 
the stick phase with increasing drive rate. As many of 
these features are common to two different systems, it is 
likely that these are general features in other stick-slip 
situations with multiple time scales that are limited to 
a window of drive rates with multiple participating time 
scales. 

A few comments may be in order about the model, 
in particular about the parameters that are crucial for 
the dynamics. While the agreement of several statistical 
and dynamical features of the model ( for several sets of 
parameter values) with the experimental AE series is en- 
couraging, it would be interesting to verify model results 
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for other sets of parameters. For instance, it is clear 
that the roller inertia and the inertia of the tape mass 
are experimentally assessable parameters. Thus, the in- 
fluence of these two inertial time scales can in principal 
be studied in experiments. However, conventional exper- 
iments have been performed keeping these parameters 
fixed, presumably, as there has been no suggestion that 
the dynamics can be sensitive to these variables. It would 
be interesting to verify the predicted dynamical changes 
in the AE signals as a function of these two parameters. 
As for the influence of the dissipation parameter 7„, the 
range of physically reasonable values of 7„ is expected 
to be small (10^'' to lO^'^) as argued. Interestingly, the 
region of low is indeed the region where both statistics 
and dynamical features compare well with that of the ex- 



periments. However, within the scope of the model, the 
visco-elastic properties of the adhesive have been mod- 
eled using an effective spring constant. (This kind of as- 
sumption is common to studies in adhesion and tackiness 
etc ^62,].) However, it is possible to include this feature as 
well. Finally, it must be stated that features that criti- 
cally depend on thickness of the film and its visco-elastic 
properties such the shape of peel front are beyond the 
scope of the present model. 
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